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DETECTION  PERFORMANCE  OF  A  SUBOPTIMUM  PROCESSOR 
IN  NON-GAUSSIAN  NOISE 

INTRODUCTION 


The  optimum  processor  for  detection  of  a  weak,  zero-mean  stationary  signal  in  stationary 
non-Gaussian  noise  was  presented  in  reference  1.  If  p{x)  is  the  first-order  probability  density 
function  (PDF)  of  the  independent  input  noise  background  samples,  and  q{x)  =  log  p{x),  then 
the  optimum  processor  for  low  input  signal-to-noise  ratio  (SNR)  takes  the  form 

K  K 

L  =  (1) 

k=\  kj=l 

where  K  is  the  number  of  samples  of  input  data  }  available  and  {p^ }  are  the  elements  of  the 

input  signal  normalized  covariance  matrix  p  (see  reference  1,  equation  (16)).  Some  of  the 

problems  associated  with  determination  of  the  q  function  and  its  first  two  derivatives  are 
discussed  in  reference  1. 

Another  problem  with  the  realization  of  processor  (1)  is  the  determination  of  the  signal 
normalized  covariance  matrix  p.  At  low  input  SNRs,  when  the  presence  of  a  signal  is  not 
known  and  the  signal  has  never  been  observed  before,  detailed  knowledge  of  matrix  p  is  very 
unlikely.  The  presence  of  signal  multipaths  of  unknown  structure,  delays,  and/or  strengths 
further  confounds  the  situation,  making  practical  realization  of  processor  (1)  virtually 
impossible.  Accordingly,  it  is  desired  to  modify  the  form  of  the  optimum  processor  to  obviate 
the  need  for  the  exact  signal  characteristics.  Also,  it  would  be  desirable  to  alleviate  the  precise 
determination  of  the  derivatives  of  the  q  function  and  replace  these  functions  with  some  simpler 
nonlinearities.  None  of  these  modifications  can  be  achieved  without  some  sacrifice  of 
detectability  performance.  Therefore,  any  suboptimum  processors  must  be  compared  with 
optimum  processor  (1)  and  a  quantitative  measure  of  losses  ascertained,  so  that  viable  decisions 
can  be  made  on  the  worthiness  of  any  replacement  candidates.  These  comparisons  must  be  done 
on  the  basis  of  the  receiver  operating  characteristics  (ROCs),  so  that  tradeoffs  in  false  alarm 
probability  P f  and  detection  probability  are  directly  quantified.  Then,  the  differences  in 

required  input  SNRs  for  equal  detection  capability  (Pj,Pj)  will  constitute  a  valid  basis  of 

comparison  of  the  various  suboptimum  processors,  and  enable  a  decision  as  to  whether  an 
acceptable  alternative  processor  exists. 
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NOISE  STATISTICS 


Three  different  noise  environments  will  be  considered — Gaussian  noise,  sech  noise,  and 
fourth-power  noise.  Descriptions  of  all  three  cases  follow. 


GAUSSIAN  NOISE 


The  relevant  statistics  of  zero-mean  Gaussian  noise  are  given  by 


p(x)  = 


X 


2  \ 


2o-' 


for  all  X, 


q{x)  =  log  p{x)  =  -log(V^  <t)  ■ 


9'w=-4,  «'w=-4. 

(7  a 


(2) 


where  cr  is  the  noise  standard  deviation.  Device  q'  is  linear  for  all  input  arguments;  the  minus 
factor  is  irrelevant  in  processor  (1).  On  the  other  hand,  nonlinearity  q”  totally  ignores  the  input 
data  {x,^}  and  can  be  dropped  from  processing  altogether,  since  additive  constants  do  not  affect 

the  ROCs  of  a  processor.  The  q  functions  of  equation  (2)  are  depicted  in  figure  1  for  <t  =  1 .  The 
optimum  processor  of  equation  (1)  takes  the  form 

K 

Z  Pv  **  (3) 

kj=l 


when  irrelevant  scale  factors  and  additive  constants  are  dropped.  Additional  properties  of 
Gaussian  noise  that  are  relevant  to  this  study  are  given  in  reference  1,  equations  (18)  through 
(27);  in  particular,  compare  equations  (19)  and  (27). 


SECH  NOISE 


Sech  noise  is  characterized  by  PDF 

PW  =  :;— sech  --  =- - - - -  for  all  ;c,  (. 

U  (t) 

where  <t  is  the  noise  standard  deviation.  This  PDF  has  exponential  decay  for  large  arguments: 

1  (  n  Uh 

^  exp  3S  X  — ^  ioo,  f  i 

<T  [  2  a) 

However,  this  PDF  is  well  behaved  at  the  origin,  in  distinction  to  exponential  noise  (see 
reference  1,  equation  (28)  and  appendix  C). 


The  corresponding  moment-generating  fonction  (MGF)  is 
=  \dx  e\p(Ax)  p(x)  = - ? -  for  \a\<—. 


PW  =  Xdx  exp{Ax)  p(x)  = -  for  U  < - ,  (f. 

^  cos(o-i)  2cr  ^ 

upon  use  of  reference  2,  equation  3.51 1  4.  Expansion  of  the  logarithm  of  the  MGF  according  to 

\ogp{A)  =  \{(7  Af  +-^(aAy  +-^((7  Ay  +  2^(crA)*  +  •••  (7 

leads  immediately  to  the  nonzero  cumulants 

2'2  =  cr',  2r4  =  2cr\  =  16c^^  Zs  =  272cr*,... .  (g 

The  c-function,  namely, 

» 

c(u,A)=  jdxexp(Ax)  p(x),  (9 

-QC 

can  be  expressed  only  in  terms  of  hypergeometric  functions. 

The  cumulative  distribution  function  (CDF)  corresponding  to  sech  PDF  (4)  is 


"f  2  r  r;r 

c(x)  =  jdu  p{u)  =  —  atan  exp - for  all  x 

'  TT  [  2  (T 
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The  inverse  of  this  CDF  follows  according  to 


2  f;r  y\ 

y  =  c(x%  x  =  c(y)  =  cr  — log  tan  —y  forO<>^<l. 

^  I  [2  J 


This  result  is  needed  to  generate  random  variables  (RVs)  with  the  PDF  in  equation 
(4)  from  uniformly  distributed  RVs  over  (0,1). 

The  q  functions  corresponding  to  PDF  (4)  are,  for  all  x. 


q(x)  =  log  p(x)  =  -  log(2  a)  -  log  cosh  —  —  , 

V  \  2  (J  j  j 


q(x)  =  -  —  -tanh  -- 
£7  2  1  2  £7 


q\x)  =  --^^  sech^l  ^  —  ]  =  qXxf 
£7  4  [2  crj 


£7^  4 


These  functions  are  depicted  in  figure  2  for  £7  =  1.  The  q'(x)  function  has  slope - - - at 

£7^  4 

1  ^ 

the  origin  and  rapidly  saturates  at  values  ± - as  |jc|  exceeds  £7;  thus,  this  nonlinearity  behaves 

£7  2 

1 

as  a  soft  limiter.  The  q^x)  function  has  a  negative  bump  at  the  origin  of  magnitude  — - - 

<7  4 

and  decays  rapidly  to  zero  for  large  x;  in  fact,  its  major  contribution  takes  place  only  for  |x| 
within  £7  of  the  origin.  There  are  no  anomalies  in  the  behavior  of  function  ^"(x)  at  the  origin, 

as  there  are  for  the  exponential  PDF  2‘‘^  exp(-2^^  |^|/cr). 


5 


The  optimum  processor  for  noise  PDF  (4)  takes  the  form 


i=l  k,M 


1 

71^ 

4 

1 

4 

where  RV 


foj. 


The  leading  scale  factor  in  the  second  and  third  lines  of  equation  (13)  can  be  dropped  without 
affecting  the  performance  of  the  optimum  processor;  however,  the  scale  factor  relating  a*  to 
must  be  retained  and  utilized.  Therefore,  the  noise  standard  deviation  cr  must  be  known. 

PDF  of  Sum  of  Independent  Seek  RVs 

The  characteristic  function  (CF)  of  a  sech  RV  is  available  from  the  MGF  in  equation  (6)  as 

p{i^)  =  ftir  exp(/  ^  x)  p{x)  = - for  all  real  (15 

*’  cos(z  cr  q)  cosh(cr  q) 

which  has  the  identical  form  as  the  PDF  in  equation  (4).  Therefore,  the  CF  of  a  sum  u  of  n 
independent  sech  RVs  with  the  same  standard  deviation  is 


cosh"(o-^)’ 

Thus,  the  PDF  of  the  sum  u  is  given  by  the  Fourier  transform 

Itt  _l  cosh"  (cr(^)  n  '  cosh"  (cr^)  na  v2  2cr  2  2cr 


2  n  lu  n  _  2 

{n-\)\7i<T  V2  2a)  v2  2crj  («-l)!;rcr  v2  2cr J 


upon  use  of  reference  2,  equations  3.512  1  and  8.380  1,  and  reference  3,  equation  6.1.23. 
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Closed  forms  can  be  developed  for  this  PDF  for  any  n.  By  using  reference  3,  equations 
6.1.15, 6.1.30,  and  6.1.31,  the  following  examples  were  developed: 

/-  \  1  1  1,  ^  ^  r- 

P\  W  =  - - 77777’  U  -  ——  (in  agreement  with  equation  (3)), 

2  a  cosh(t/)  2  cr 

,  ,  1  U 

P2(«)  = - rT7777’ 

Tta  smh(t/) 


4cri^  j  cosh(t/) 

/  ^  1  r.  ^  f/ 

P,i.u)  =  - -  4  +  —  — — , 

OTta  y  a  j  smh(t/) 

,,  1  fo  1 

48(7  a^  jy  cr^  J  cosh(t/) 

m;ra[  a^ )[  (j^  I  sinh(t/)' 


In  fact,  a  simple  recursive  formula  for  the  wth  PDF  follows  readily  from  equation  (17): 

.  .  (n-iy  +u^/a^  . 

Pn(M)  =  P„-2(m)— — 777 - ^  for«>3. 

Starting  values  are  available  from  equation  (18). 


FOURTH-POWER  NOISE 

The  PDF  for  this  type  of  noise  is  given  by 


Pix)  = 


3/(V2;r) 


1  +  —^  1  + - - 

I  a  II  2cr^ 


a^lTT  ^  ^  2x^ 


for  allx. 


where  a  is  the  noise  standard  deviation.  For  large  arguments,  this  PDF  decays  only  as  1  /x\ 
which  constitutes  a  very  heavy-tailed  (HT)  distribution.  The  corresponding  CF  is 

/^0'^)  =  2  exp(-|^|<T/V2)-exp(-|4^|(TV2)  forallreal^.  (; 


The  c(m.  A)  function  of  equation  (9)  involves  four  exponential  integrals. 
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The  CDF  for  fourth-power  PDF  (20)  is 


11  f  R  \  f  \  (  (  f 

,,11-  ^J2  X  X  2 

c(x)  =  —  +  —  2  atan  -  -  atan  —j= —  =  —  atan  exp  3  asinh  —j= — 

2  71  I  cr  J  I V2  cr  I  7t  i  V2  cr 


The  inverse  to  this  relation,  y  =  c(jc),  is  given  by 


|1  f  ( 71  HI 

;c  =  c(>’)  =  cr-\/2sinh  —log  tan  —y  for0<y<l. 
13  I  V  2  J) 


(22) 


(23) 


The  derivation  of  the  last  two  relations  is  given  in  appendix  A.  Result  (23)  allows  for  rapid 
generation  of  RVs  with  fourth-power  PDF  (19)  directly  from  uniformly  distributed  RVs  over 
(0,1). 

The  q  functions  for  fourth-power  PDF  (20)  are 


These  functions  are  depicted  in  figure  3  for  cr  =  1.  Function  q’  is  linear  for  small  arguments, 
but  then  saturates  and  returns  slowly  to  zero  as  -Alx  for  large  x.  The  second  derivative  q"  has  a 
large  negative  bump  near  the  origin  and  returns  to  zero  as  4/x^  for  large  jc.  The  optimum 
processor  is  available  upon  reference  to  equation  (1). 

Plots  of  the  three  types  of  noise  processes  are  given  in  figure  4  for  a  standard  deviation  of 
rr  =  1 .  Whereas  the  Gaussian  sample  data  are  limited  to  a  range  of  ±3,  the  sech  noise  sample 
reaches  values  of  ±5,  and  the  fourth-power  noise  sample  attains  values  near  ±10.  Also,  the  burst¬ 
like  character  of  the  fourth-power  noise  is  very  obvious,  occasionally  reaching  values  many 
times  that  of  its  standard  deviation.  It  is  these  bursts  that  cause  serious  degradation  of 
performance  of  processors  designed  to  operate  in  a  Gaussian  noise  background. 
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q  FUNCTIONS  FOR  NONLINEARLY  TRANSFORMED  RANDOM  VARIABLES 

The  three  noise  processes  discussed  above  were  generated  directly  from  uniformly 
distributed  RVs  over  (0,1).  An  alternative  method  of  generating  RVs  is  by  means  of  nonlinear 
transformations  of  RVs  with  already  established  statistics  (other  than  uniform).  Here,  RV  y  will 
be  generated  from  RV  x  according  to  monotonically  increasing  transformation yi  y  =y(x). 
Random  variable  x  has  PDF  p^{x)  and  corresponding  (zeroth-order)  q  function 

q^{x)  =  \ogp^{x).  (25: 

The  first  nonlinear  transformation  of  interest  is  cubic,  resulting  in  RV 

y  =  x  +  flx\  cr>0  (for monotonicity).  (261 


y  =  f{x)  =  x  +  ax\  /'(x)  =  l  +  3^^x^  x^f(y)  =  F, 
where  f  denotes  the  inverse  function  to  f.  This  inverse  function  is  given  by 


P  =  f{y)^t{y)--  I  t{y)= 

2>at{y)  la  vUaJ  21 


(see  reference  1,  equations  (42)  through  (46)).  Then,  upon  use  of  equation  (27), 


F'  =  f\y)  = 


/'(/(T))  nn  l  +  3aF^ 


This  result  affords  a  simple  technique  of  finding  higher-order  derivatives  of  inverse  function  F. 
The  PDF  of  RV  y  is  given  by 

p.(f) 

/'(/w)  nn  i+3aF=' 

As  an  example  of  an  application,  let  RV  x  be  zero-mean  Gaussian  with  standard  deviation 
a .  Then,  the  PDF  of  y  =  x  +  a  xMs 


pAy)  = 


1  exp(-FV(2rT^)) 

■yjln  <7  \  +  2>aF^ 


The  first  three  q  functions  follow  immediately  upon  use  of  equations  (29)  and  (30)  as 
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qy  (j)  =  log  (y)  =  -  logiy/l^  a)  -  -  log(l  +  3  a  ), 


q'yiy)  =  - 
q^yiy)  =  - 


\  F  ^  F 
cr'  l  +  3aF^  {\  +  3aF^f 

1  \-3aF^  \-9aF^ 

cr'  (l  +  3aF')'  ^  {\  +  3aFy 


(32) 


These  equations  quickly  confirm  the  numerical  results  in  reference  1,  figure  3,  which  were 
obtained  by  a  much  more  laborious  technique  involving  symbolic  differentiation.  The  function  F 
is  given  by  equation  (28)  above.  These  results  allow  for  realization  of  the  optimum  processor  in 
equation  (1). 

Plots  of  the  q  functions  in  equation  (32)  are  displayed  in  figures  5,  6,  7,  and  8  for  scale 
factor  a  equal  to  0.1,  0.2,  0.3,  and  0.4,  respectively  (a  =  0  corresponds  to  a  Gaussian  RV).  The 
four  g  functions  in  the  uppermost  subplots  gradually  develop  a  more  pronounced  flare  for  large 
arguments as  a  increases.  Also,  the  q”  functions  develop  a  progressively  larger  negative 
bump  near  the  origin.  The  gradual  narrowing  of  the  various  functions  in  y  in  the  figures  is  due  to 
the  standard  deviation  of  RV  y  increasing  as  a  increases;  thus,  the  standard  deviations  of  y  for 
the  scale  factors  detailed  above  are  1.32, 1.67, 2.04,  and  2.41,  respectively.  That  is,  abscissa 
value y  =  1  in  the  figures  does  not  correspond  to  unit  standard  deviation,  as  it  did  in  the  earlier 
figures. 

As  j  ->  +00,  the  asymptotic  behaviors  of  the  functions  above  are  according  to 


'O’) -(Wo)"’,  F-(ylaf\  J 


q'yiy) 


-1 


1 


3«''V 

/?y(y)  =  exp(^y(>;))~exp| 


qliy)' 


2<y^ 

2/3  2  ..4/3 


9a"' V"  y 


(33) 


1  1 

f-1 

2/3\ 

2cr^l 

UJ 

) 

Thus,  both  derivatives  of  the  q^iy)  function  decay  to  zero  for  large  arguments,  albeit  at  very 

slow  rates,  especially  for  the  first  derivative.  The  decay  of  the  PDF  of  y  is  slower  than 
exponential.  All  the  results  above  apply  only  to  a  Gaussian  RV  x  and  cubic  nonlinear 
transformation  (26). 
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More  generally,  for  an  arbitrary  monotonically  increasing  nonlinear  transformation/ and 
input  PDF  p^,  the  relevant  transformations  are  defined  according  to 

y  =  f(x),  X  =  g(y),  ^4) 

where  g  is  the  inverse  function  to/  It  is  presumed  that  function /and  its  derivatives  can  be 
readily  evaluated.  The  CDF  of  RV  y  =  /(x)  is  given  by 

S  (t)  =  Pr(y  <y)  =  Pr(/(x)  <y)  =  Pr(x  <  ^(y))  =  c,  (g(y)).  (35) 

The  corresponding  PDF  of  y  is 

Pyiy)  =  c;iy)  =  c[{giy))  g'(y)  =  pMy))  gXy)-  (36) 

But,  from  equation  (34), 


y  =  figiy)),  \  =  f'{g(y))g'iy),  g'(y)  =  — (37 

/  (g(y)) 

This  relation  avoids  the  need  to  evaluate  the  derivative  of  the  inverse  function  g.  Substitution  of 
equation  (37)  in  equation  (36)  yields 


Pyiy) 


p.jgjy)) 
f'igiy))  ■ 


(38) 


The  zeroth-order  q  function  for  RV  y  is 

(iy  (T)  =  log  Py  iy)  =  log  p,  (g(y ))  -  log  f’ig(y))  =  (g(y))  -  log  f'ig{y)).  (39) 

The  first  derivative  follows  as 


_  (I'Agjy))  f\giy)) 
f'igiy))  f'igiy)?' 


(40) 


where  equation  (37)  has  been  used.  The  second  derivative  follows  in  a  similar  manner  as 
n”(v\  =  <^Siy))f'igiy))-q:igiy))  figiy))  rigiy))  f\g{y))  -  2  rigiy)?  , 

The  functions  needed  are  the  q  functions  for  RV  x,  the  transformation/and  its  first  three 
derivatives,  and  the  inverse  function  g. 

If  inverse  function  g  cannot  be  determined  analytically,  the  following  numerical  approach 
will  yield  sample  values  of  the  desired  q  functions  of  RV  y.  Select  a  set  of  fine  spacings  in  x. 
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namely,  }  for  «  =  1 :  iV,  where  N  is  large  enough  to  cover  the  range  of  interest.  Then,  from 
transformation  (34),  evaluate  and  store  the  values 

=  fix„ )  forn  =  l:N;  equivalently, =  g(y„ ).  (42) 

Now,  equations  (39)  through  (41)  immediately  yield,  for  n  =  l:  N, 


9'y(T«)  =  9x(jf„)-log/'(JC„), 

’  "  f'M" 

’  fix.?  fXxj' 


All  of  the  quantities  needed  here  are  readily  available  from  the  tabulation  in  equation  (42). 
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PERFORMANCE  OF  OPTIMUM  PROCESSORS  IN  NON-GAUSSIAN  NOISE 


The  background  noise  will  be  taken  to  be  one  of  the  three  types  introduced  in  the  previous 
section,  namely,  Gaussian  noise,  sech  noise,  or  fourth-power  noise.  The  signal,  when  present, 
will  be  the  multipath  model  given  in  reference  1,  page  9.  Namely,  signal  sample 

s*  =  ^  ^k-m  ’  w  >  0,  (44) 

for  k  =  \  -.K,  where  zero-mean  sequence  {d  J  consists  of  uncorrelated  samples.  In  all  the 
following  numerical  results,  the  number  of  independent  noise  data  samples  \s.  K  -  1000,  the 
signal  delay  is  w  =  10,  and  the  relative  attenuation  factor  is  or  =  -1,  corresponding  to  a  perfectly 
reflecting  surface-bounce  path.  The  ROCs  have  been  generated  by  using  one  million  trials  for 
each  processor  and  signal-noise  combination  of  interest.  The  input  SNR  is  the  ratio  of  the  total 
power  in  signal  sample  (44)  to  the  noise  power  per  sample. 

The  initial  ROCs  in  figure  9  pertain  to  the  usual  reference  case,  namely,  the  optimum 
processor  for  Gaussian  signal  in  Gaussian  noise  (see  equation  (3)  and  reference  1,  equation  (35)). 
Equation  (3)  will  be  referred  to  as  Processor  G.  Signal  parameters  a  and  m  are  presumed 
known  to  the  receiver  processor.  Three  operating  points  (OPs)  are  of  particular  interest:  low- 
quality  OP  Pj  =  le-3,  Pj  =  0.5;  medium-quality  OP  Pf  =  le-4,  P^  =  0.9;  and  high-quality  OP 

Pf  ^U-5,Pj  =  0.99.  The  required  input  SNRs  to  realize  these  three  OPs  are  observed  to  be 

-9.2  dB,  -6.8  dB,  and  -5.3  dB,  respectively.  These  numerical  results  are  consistent  with  the 
original  derivation  in  reference  1,  which  was  for  low  input  SNRs. 

The  ROCs  in  figure  10  differ  only  in  that  the  background  is  changed  to  sech  noise,  while  the 
processor  is  kept  as  equation  (3),  with  all  signal  parameters  known,  namely,  Processor  G.  This  is 
not  the  optimum  processor  for  a  sech  noise  background.  The  required  SNRs  at  the  low-, 
medium-,  and  high-quality  OPs  are  —7.85  dB,  —5.6  dB,  and  —4.1  dB,  respectively.  These  values 
are  not  significantly  different  from  those  in  figure  9,  being  approximately  1 .2  dB  worse. 

However,  when  the  background  is  changed  to  fourth-power  noise  and  the  processor  is  still 
unchanged  from  that  above,  the  ROCs  indicate  grossly  degraded  performance,  as  shown  in  figure 
11.  The  required  input  SNR  to  realize  the  low-quality  OP  is  approximately  +6  dB,  a  whopping 
15-dB  degradation  relative  to  the  Gaussian  background  case.  This  example  demonstrates  the 
dangers  in  trying  to  use  the  processor  designed  for  Gaussian  interference,  namely.  Processor  G, 
in  the  presence  of  HT  noise,  even  when  the  signal  parameters  are  exactly  known.  The  losses  for 
this  suboptimum  processor  at  the  medium-  and  high-quality  OPs  are  still  worse  and  were  not 
estimated. 


To  determine  bases  of  performance  in  non-Gaussian  noise  backgrounds,  the  optimum 
processors  for  those  environments  must  be  investigated  quantitatively.  The  optimum  processor 
for  a  sech  noise  background  was  presented  in  equation  (13)  and  is  denoted  as  Processor  S.  The 
optimum  processor  for  a  fourth-power  noise  background  is  available  upon  combination  of 
equations  (1)  and  (24)  and  is  denoted  as  Processor  F.  These  optimum  processors  require  and 
utilize  knowledge  of  the  signal  characteristics  and  parameters. 
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To  determine  the  fundamental  capability  of  detection  in  sech  noise,  and  to  ascertain  the 
dependence  on  the  particular  signal  type,  the  ROCs  in  figures  12  through  14  were  generated.  All 
three  cases  utilize  Processor  S,  while  the  signal  is  alternatively  Gaussian,  sech,  and  fourth-power, 
respectively.  The  required  input  SNRs  are  rather  similar  for  the  Gaussian  and  sech  signal  cases, ' 
while  those  for  the  fourth-power  signal  are  about  1 . 1  to  1 .7  dB  larger.  Thus,  the  optimum 
detectability  does  vary  slightly  with  the  signal  type,  with  the  Gaussian  signal  being  most 
detectable  and  the  fourth-power  signal  being  the  least  detectable.  The  HT  character  of  the 
fourth-power  signal  does  not  help  detectability  for  low  input  SNRs. 

The  corresponding  results  for  optimum  detection  in  fourth-power  noise,  that  is,  for  processor 
F,  are  presented  in  figures  15  through  17.  The  three  sets  of  ROCs  correspond  to  Gaussian,  sech, 
and  fourth-power  signals,  respectively.  Again,  the  required  input  SNRs  for  the  Gaussian  and 
sech  signal  cases  are  rather  similar,  while  those  for  the  fourth-power  signal  are  about  1 .2  to 
2.2  dB  larger,  depending  on  the  exact  operating  point  of  interest.  Alternative  operating  points 
can  easily  be  investigated  from  the  ROCs  presented  in  figures  9  through  17. 
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DERIVATION  OF  SUBOPTIMUM  PROCESSORS 


The  form  of  the  optimum  processor  was  presented  in  equation  (1)  and  is  repeated  here: 

K  K 

L  =  £^''(Xi)+  (45) 

*=1  kj=\ 

Its  realization  requires  knowledge  of  nonlinearities  q'  and  q”  and  the  signal  normalized 
covariance  matrix  p.  An  alternative  expression  follows  as 


K  K  K 

L  =  Z9''(Xi)  +  X^'(*t)^  ^'(Xi)9'(x,)-  (46) 

*=1  /t=l  k,j=\ 

l‘*j 


All  the  off-diagonal  terms  in  equation  (46)  require  information  about  the  signal  characteristics. 
In  the  absence  of  this  signal  knowledge,  dropping  these  terms  suggests  a  much  simpler  but 
suboptimum  processor,  namely. 


(47) 


The  degradation  of  performance  caused  by  this  significant  reduction  depends  on  the  detailed 
signal  characteristics  and  the  particular  background  noise  statistics. 

The  difficulty  of  accurately  estimating  the  q'  and  q”  functions  fi'om  limited  amounts  of 
noise-alone  data  also  begs  for  some  approximation  of  the  nonlinearities  indicated  in  equation 
(47).  A  clue  to  a  possible  simplification  is  afforded  by  the  sech  noise  results  in  the  third  line  of 
equation  (12).  Namely,  is  the  square  of  the  additive  constant  can  be  dropped  since  it 
does  not  affect  the  ROCs.  Furthermore,  the  second  line  of  equation  (12),  along  with  figure  2, 
indicates  that  is  a  soft  limiter.  That  is,  the  suboptimum  processor  in  equation  (47)  can 
reasonably  be  replaced  by  a  soft  limiter  followed  by  an  energy  detector: 


^2 

ji=i 


yk 


f 

=  ±tanh 

V 


^  A 
—  ** 

2  fj  ; 


for  A:  =  \:K. 


(48) 


(The  polarity  of  is  immaterial,  because  it  is  squared  before  summation.)  Notice  that 
knowledge  of  the  noise  standard  deviation  <j  is  needed  for  realization  of  processor  (48). 

More  generally,  inspection  of  figures  2,  3,  and  5  through  8  reveals  that  the  general  behavior 
of  nonlinearity  q'  is  that  of  a  soft  limiter,  while  that  of  q”  is  similar  to  the  square  of  q',  less  an 
additive  constant.  This  suggests  still  a  further  simplification  to  the  suboptimum  processor, 
namely,  a  hard  limiter  followed  by  energy  detection: 
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for  =  1 :  a:. 


(49) 


i=l 


y*  = 


[6  sgn(x J  for|xJ>6j 


The  parameter  b  must  be  chosen  proportional  to  the  noise  standard  deviation  a  to  realize  the 
best  detectability  performance.  (This  operation  is  distinct  from  clipping,  for  which  sgn(X;t ) 
would  be  used  for  all  arguments.)  The  overall  effect  of  the  hard  limiting  and  squaring  is 
equivalent  to  subjecting  the  input  data  samples  {x^^ }  to  a  squaring  operation  and  saturation: 


K 

k=l 


A  for|xJ<6 
for|xJ>A 


for  A:  =  1 :  /C. 


(50) 


Finally,  since  additive  constants  do  not  affect  the  ROCs,  this  processor  can  be  modified  to 


K 


\A-b^  for|x;t|<6| 
0  for|x;t|>6 


k=\ 


=  min(x^  -  b^  ,0)  for  A  =  1 : 


(51) 


This  nonlinearity  applied  to  the  input  data  is  simply  a  negative  parabolic  bump  centered  about 
the  origin,  and  zero  for  input  magnitudes  greater  than  b.  Reconsideration  of  figures  2, 3,  and  5 
through  8  reveals  that  is  a  very  reasonable  approximation  to  the  fimction  q".  Furthermore, 

processor  L4  has  been  deduced  through  an  argument  that  employed  a  plausible  and  physically 
reasonable  processor  in  HT  noise,  namely,  hard  limiting  followed  by  energy  detection.  This 
argument  also  lends  physical  interpretation  to  the  q"  term  in  optimum  processor  (1).  The  only 
noise  parameter  that  must  be  monitored  is  the  noise  standard  deviation;  this  is  a  much  simpler 
and  stable  task  than  estimating  the  entire  second  derivative  of  the  noise  PDF.  How  well  this 
suboptimum  processor  performs  in  HT  noise  distributions  must  now  be  determined  by 
simulations. 
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PERFORMANCE  OF  A  SUBOPTIMUM  PROCESSOR 


The  first  case  of  interest  is  suboptimum  processor  (51)  for  a  sech  noise  background.  The 
ratio  of  saturation  parameter  b  to  the  noise  standard  deviation  g  is  denoted  by  ratio 

b 


Through  trial  and  error,  the  best  value  of  ratio  r,  when  the  background  is  sech  noise,  was  found 
to  be  approximately  1,  regardless  of  the  signal  type.  This  value  of  r  maximized  the  detection 
probability  for  a  given  false  alarm  probability  Pj.  The  signal  model  is  the  same  as  that 

described  in  equation  (44)  and  the  accompanying  discussion.  The  ROCs  for  Gaussian,  sech,  and 
fourth-power  signals  in  sech  noise  are  presented  in  figures  1 8  through  20,  respectively.  The 
losses  of  this  suboptimum  processor  relative  to  the  optimum  results  in  figures  12  through  14  are 
given  in  table  1 . 


Table  1.  Losses  of  Suboptimum  Processor  in  Sech  Noise 


Signal  Type 

Losses  (dB) 

Gaussian 

1.0 

0.95 

0.95 

Sech 

1.0 

0.95 

0.85 

Fourth  Power 

1.05 

1.0 

1.0 

These  are  very  slight  losses  in  detectability,  considering  the  facts  that  the  signal  structure  is 
ignored  by  suboptimum  processor  (51)  and  that  the  noise  nonlinearities  q'  and  q”  were  not 
estimated  or  employed.  In  fact,  the  uniformity  of  the  losses  over  the  range  of  signal  statistics 
(three  very  different  types)  and  the  range  of  operating  points  (three  widely  differing  qualities)  is 
surprising. 


V 
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The  corresponding  ROCs  for  fourth-power  noise  are  presented  in  figures  21  through  23. 
Here,  the  best  value  of  ratio  r  in  equation  (52)  was  found  to  be  approximately  0.65,  regardless  of 
the  signal  type.  The  losses  of  the  suboptimum  processor  (51)  relative  to  the  optimum  results  in 
figures  15  through  17  are  given  in  table  2. 

Table  2.  Losses  of  Suboptimum  Processor  in  Fourth-Power  Noise 


Signal  Type 

Losses  (dS 

0 

Gaussian 

0.8 

0.8 

0.8 

Sech 

0.9 

0.85 

0.7 

Fourth  Power 

0.8 

0.7 

0.5 

The  losses  for  fourth-power  noise  are  less  than  those  for  sech  noise  in  table  1.  However,  the 
results  in  table  2  are  slightly  more  variable  with  regard  to  signal  type  and  quality  of  the  operating 
point.  Nevertheless,  the  overwhelming  feature  of  tables  1  and  2  is  the  relatively  insignificant 
losses  associated  with  employment  of  suboptimum  processor  (51).  Furthermore,  this  comparison 
has  been  made  with  the  absolute  optimum  processor  (1),  which  requires  knowledge  of  the 
detailed  signal  structure;  this  signal  knowledge  will  usually  not  be  known  in  sufficient  detail  to 
make  much  use  of  it  in  practice. 

Several  additional  values  of  ratio  r  in  equation  (52)  were  simulated  in  addition  to  the  value 
0.65  cited  above  for  fourth-power  noise.  The  ROCs  were  not  very  sensitive  to  the  exact  value  of 
r  employed.  An  example  for  r  =  0.75  is  given  in  figure  24;  the  performance  is  only  slightly 
poorer  than  in  figure  21.  However,  the  excessive  simulation  time  required  to  try  the  various  r 
values  prompted  the  consideration  of  a  simpler  detectability  measure.  In  appendix  B,  the 
deflection  criterion  is  investigated  for  Gaussian,  sech,  and  fourth-power  noise  and  suboptimum 
processor  (5 1).  The  numerical  results  on  the  deflection  criterion  confirm  the  choices  for  r  made 
above  with  respect  to  the  ROCs  themselves.  This  alternative  approach  in  appendix  B  can  be 
used  to  set  the  ratio  r  “on  the  fly”  in  a  practical  setting  where  lengthy  ROCs  could  not  be 
computed. 
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0.999 

0.990 


40 


Thus  far,  the  signal  model  has  been  that  of  equation  (44)  with  relative  attenuation  factor 
a  = -  I  and  m  =  10.  To  ascertain  the  dependence  on  the  signal  structure,  the  factor  a  is  now  set 
equal  to  —0.5.  The  received  signal  is  Gaussian  and  the  background  is  maintained  as  fourth-power 
noise.  The  ROCs  for  the  optimum  processor  are  presented  in  figure  25,  while  those  for  the 
suboptimum  processor  are  displayed  in  figure  26.  When  factor  a  is  set  to  zero  (no  surface 
reflection,  making  m  irrelevant),  the  corresponding  ROCs  for  the  optimum  and  suboptimum 
processors  are  given  in  figures  27  and  28,  respectively. 

The  dependencies  of  the  detectability  losses  of  the  suboptimum  processor  relative  to  the 
optimum  processor  are  listed  in  table  3.  The  top  line  comes  from  table  2,  while  the  bottom  two 
lines  come  from  figures  25  through  28.  The  three  numbers  on  a  line  correspond  to  the  low-, 
medium- ,  and  high-quality  OPs,  respectively.  Two  features  stand  out:  the  variation  with  OP  is 
virtually  nil,  and  the  degradation  with  respect  to  relative  attenuation  factor  a  decreases  to  zero 
as  a  does.  This  is  not  precisely  true;  however,  there  is  very  little  difference  between  the  ROCs 
in  figures  27  and  28. 


Table  3.  Dependence  of  Losses  on  Signal  Structure 


a 

Losses  (dB) 

-1.0 

0.8 

0.8 

0.8 

-0.5 

0.6 

0.5 

0.4 

0.0 

0 

0 

0 

Finally,  a  signal  with  more  structure,  rather  than  less,  was  considered.  Signal  model  (44) 
was  expanded  to  include  three  paths  with  relative  strengths  +1,-1,  +1  and  delays  0,  10,  20, 
respectively.  The  ROCs  for  the  corresponding  optimum  processor  are  given  in  figure  29,  while 
the  ROCs  for  the  suboptimum  processor  are  given  in  figure  30.  Now  the  loss  is  a  uniform  1 .5  dB 
across  all  the  OPs.  This  loss  is  greater  than  any  of  those  encountered  in  tables  1  through  3  and  is 
due  to  the  increased  off-diagonal  structure  of  the  signal  covariance  matrix.  This  increasing  loss 
can  be  expected  to  prevail  as  the  off-diagonal  structure  increases  in  complexity.  However,  the 
likelihood  of  knowing  this  information  in  advance  is  not  reasonable  in  practice.  Therefore,  the 
relatively  small  losses  listed  above  must  be  taken  as  extremely  encouraging  in  detection  of  a 
weak  signal  in  HT  noise. 
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Figure  26.  R 
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SUMMARY 


The  optimum  processor  has  been  derived  for  optimum  detection  of  a  weak  stationary  signal 
with  multipath  in  the  presence  of  additive  non-Gaussian  noise.  The  optimum  processor  requires 
knowledge  of  the  signal  normalized  covariance  matrix  and  the  first  two  derivatives  of  the  first- 
order  probability  density  function  of  the  noise.  The  difficulty  in  obtaining  some  or  all  of  this 
information  in  practice  suggests  alternative  simpler  processors  that  require  only  partial  or 
minimal  information  for  their  realization.  However,  these  alternative  processors  have  degraded 
detection  capability  relative  to  the  optimum  processor,  and  the  exact  amounts  of  these  losses  are 
important  and  must  be  determined  before  any  alternatives  can  be  considered  acceptable. 

A  particular  suboptimum  processor  has  been  investigated  here,  namely,  a  hard  limiter 
followed  by  an  energy  detector.  The  hard  limiter  prevents  the  heavy-tailed  noise  contributions 
from  reaching  the  processor  output,  while  the  long-term  energy  detection  allows  for  the  weak 
signal  (if  present)  to  contribute  to  the  device  output,  enabling  its  detection  even  in  low  input 
signal-to-noise  ratio  situations.  Signal  structure  information  is  ignored  by  this  processor.  The 
losses  in  performance  typically  run  on  the  order  of  1  dB,  regardless  of  whether  low-quality  or 
high-quality  detection  performance  is  required.  This  rather  low  loss  is  very  encouraging, 
especially  in  the  presence  of  heavy-tailed  noise  distributions,  which  can  lead  to  catastrophic 
losses  in  the  usual  linear  processor  designed  for  a  Gaussian  background.  The  plausibility  of  this 
type  of  processor  is  also  confirmed  by  deriving  it  as  an  approximation  to  the  optimum  processor. 

This  quantitative  study  provides  a  basis  for  making  some  conclusions  about  operating  in  a 
non-Gaussian  environment;  such  studies  are  few  and  far  between.  It  has  concentrated  on 
accurate  evaluation  of  the  receiver  operating  characteristics  of  the  various  processors,  thereby 
furnishing  hard  numerical  signal-to-noise  ratio  comparisons.  The  cost  of  ignoring  pertinent 
information  or  using  inaccurate  estimates  can  be  accurately  assessed,  and  fiirther  studies  can  be 
initiated  to  gather  the  pertinent  required  data,  with  full  knowledge  of  the  gains  to  be  accrued. 
Also,  the  ultimate  detectability  of  a  signal  in  that  particular  non-Gaussian  environment  can  be 
ascertained,  once  and  for  all,  for  possible  future  applications. 
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APPENDIX  A 

INVERSION  OF  CDF  FOR  FOURTH-POWER  RANDOM  VARIABLE 


The  essential  problem  in  inverting  equation  (22)  in  the  main  body  of  this  report  for  the 
fourth-power  cumulative  distribution  function  (CDF)  is  the  solution  of  the  following  equation  for 
X  in  terms  of^^: 


y  =  2  atan(;c)  -atan(jc/2). 


(A-1) 


The  appropriate  steps  follow: 


y  +  atan(xl 2)  =  2  atan(Af),  tan(  y  +  atan(x/2))  =  tan(2  atan(x)). 


(A-2) 


Let  z  =  tan(>^).  Then,  equation  (A-2)  yields 


z  +  x/2  2x 
l-zx/2  1-x 


(2z  +  x)(l- x^)  =  2x(2- z x),  x^  +  'ix-2z  =  0. 


(A-3) 


Let  x  =  t  +  u  to  get 


r  +u^  -  2z  +  3(t  +  u){l  +  tu)  =  0. 


(A-4) 


Now  let  u  =  -\/t  to  get 


t^--^-2z  =  0,  t^^2zt^-l  =  0,  t^=z  +  ^J\  +  z\ 

t  =  {z  +  yl\  +  z^^  =  exp  —  log(z  +  Vl  +  )  =  exp  —  asinh(z)  . 

\3  y  V  3  J 


(A-5) 


Therefore,  the  real  solution  to  equation  (A-1)  is 


x  =  t-^,  /  =  (z-i- +  )  ,  z  =  tan(>^). 


(A-6) 


Alternative  forms  fort  directly  in  terms  of  y  are 


1  +  sin(>’) 
cos(y) 


J  ^  y 

=  cot - — 

U  2 


(A-7) 


A-1 


The  nonlinear  transformation  that  yields  zero-mean,  unit-variance  fourth-power  random 
variables  (RVs)  x  directly  from  uniformly  distributed  RVs  u  over  (0,1)  is  given  by  equation  (23) 


X  =  V2  sinh  -  log  tanf — u 
3  12 


(A-8) 


The  CDFs  and  probability  density  functions  (PDFs)  of  the  unit-variance  components  of  this 
transformation  are 

^  2  2/ 

X,  =  tan  u  ,  Cj  (x)  =  atan(A:),  (x)  = - -  for  0  <  x  <  00; 

J  K  1+X- 

2  2  {( n  \\  1  ) 

X,  =  — log(x,),  C2(jr)  =  — atan  exp  —  X  ,  /),(x)  =  --sech  —  a:  forallA:; 

^  ^  I  l2  jJ  ^  2  i2 


r  .  fn  \  2  f  f  xW) 

x=v2sinh  —  ,  c(a:)  =  — atan  exp  3asinh  -= 

J  ^  [  lV2jJJ 

for  all  X. 

(l  +  2x')(l  +  xV2) 


(A-9) 


All  of  the  results  above  and  in  equations  (20)  through  (24)  are  special  cases  of  the  following 
equations,  when  b  is  set  equal  to  3.  The  general  CDF  is 


2  X  xM  2  (  f  fxW\ 

c(x)  =  — atan  -  +  ,1  +  ^  =  -  atan  exp  6  asinh  -  for  all  x, 

^  U  V  flM  ;r  [a 


(A- 10) 


where  b  need  not  be  integer.  Letting  y  =  c(x),  the  inverse  relation  is 


ri  r  /^;z- 

x  =  c(y)  =  a  sinh  -log  tan  —  for0<j<l, 

\b  V  v2  Jl) 


(A-11) 


The  general  PDF  is  the  derivative  of  equation  (A- 10);  namely. 


.  ,  2  bj  a 

p(x)  = 


n  q  {q  +  xia)  -^(^-x/a)' 


>  9  =  +  fl>0,  6  >  0. 


(A-12) 


The  asymptotic  behavior  of  the  PDF  is  according  to  the  inverse  b+\  power: 


.  .  A  b  \  a  \ 

P\X)  — j  r  as  X  — ^  ioo. 

n  a  2\x\  i 


(A-13) 


A-2 


Thus,  RVs  with  any  desired  asymptotic  power- law  behavior  can  be  generated  by  means  of 
equation  (A-1 1),  wheny  is  replaced  by  a  uniformly  distributed  RV  over  (0,1). 

For  the  variance  to  be  finite,  that  is, 

=  ^dx  p{x)  <  00,  (A- 14) 

the  parameter  b  must  be  greater  than  2.  Then, 

cr^^a^Qib),  (A- 15) 

where 


b>2. 


(A-16) 


For  a  given  value  of  b,  relations  (A- 15)  and  (A-16)  give  the  required  value  of  parameter  a; 
namely. 


(A- 17) 


in  terms  of  the  desired  standard  deviation  cr.  The  integral  in  equation  (A-16)  can  be  evaluated  in 
closed  form  for  b  equal  to  small  integers: 

0(3)  =  j,  0(4)  =  ^^,  0(5)  =  0(6)  =  ^^^.  (A-18) 

2  2  2  6 

Additional  non-integer  b  values  will  require  numerical  integration  of  equation  (A-16).  However, 
one  special  case  that  can  be  evaluated  in  closed  form  is  0(2.5)  =  -\/5/2. 


A-3  (A-4  blank) 


APPENDIX  B 

DEFLECTION  CRITERION  FOR  SUBOPTIMUM  PROCESSOR 


The  decision  variable  is  given  by  a  sum  of  K  independent  samples  of  noise  plus  signal  (if 
present)  subjected  to  an  arbitrary  nonlinear  transformation  f: 


w  =  X/(x^)  =  Xy(A:). 


(B-1) 


The  relevant  deflection  criteria  are  defined  as 


,  d„=^d=y[K 


E(y,^„)-E(yJ 

{E(>l)-E(y,yy“’ 


where  E  denotes  an  ensemble  average.  The  averages  needed  are  given  by 


E(y„)=  jdu  f(u)  p„(u), 
E(yl)=  jdufiuY  p„{u), 
E{y,^„)=  \du  f{u)  p^„{u). 


(B-2) 


(B-3) 


To  determine  the  probability  density  function  (PDF)  p^„  when  signal  and  noise  are  present, 

a  simplifying  candidate  signal  PDF  will  be  employed;  this  is  not  expected  to  be  too  critical  for 
large  K,  when  low  input  signal-to-noise  ratios  (SNRs)  are  relevant.  In  particular,  the  signal  PDF 
is  taken  as  flat  over  interval  (-v,v),  where  There  follows 

I  V  2 

Ps. («)  =  Ps (x)  P„(u-x)  =  —  \dx  p„iu-x)  =  —  [c„ (m  +  v) -  c„ (u  -  v)],  (B-4) 

^  7  V  J  7  \/ 


where  c„  is  the  noise  CDF.  Then,  for  use  in  equation  (B-2),  along  with  equation  (B-3), 


^(y /(m)  (m  +  v)  -  c„  (m  -  v)l,  V  =  V3  (T, , 


(B-5) 


which  can  be  numerically  evaluated  for  the  nonlinearity / of  interest.  Figure  B-1  presents  a 
program  for  this  computation.  The  three  noise  PDFs  are  again  Gaussian,  sech,  and  fourth  power. 
The  nonlinearity  of  interest  is  the  suboptimum  processor  given  in  equation  (51)  in  the  main 
report,  namely. 


/ (m)  =  min(w^  -  6^  ,0). 


(B-6) 


B-1 


The  deflection  of  output  w  of  processor  (B-1 )  for  a  Gaussian  noise  background  is 

displayed  in  figure  B-2.  The  abscissa  r  is  the  ratio  b!  cr  of  the  saturation  parameter  to  the  noise 
standard  deviation.  As  r  increases,  the  deflection  increases  monotonically.  This  is  consistent 
with  the  fact  that  the  optimum  processor  in  Gaussian  noise  is  to  let  b  tend  to  infinity,  that  is,  no 
saturation  at  all.  Alternatively,  when  interpreted  as  a  hard  limiter  followed  by  an  energy 
detector,  the  hard-limiting  saturation  point  is  allowed  to  tend  to  infinity,  leading  to  the  standard 
energy  detector  for  a  Gaussian  noise  background. 

The  deflection  criterion  for  sech  noise  is  given  in  figure  B-3.  For  low  input  SNR,  the  best 
value  of  ratio  r  is  1 ;  this  is  the  value  used  for  the  ROCs  in  figures  1 8  through  20  in  the  main 
report.  This  best  value  was  determined  for  the  ROCs  by  simulation  for  several  different  values 
of  r;  here,  it  has  been  deduced  from  the  simpler  deflection  criterion.  However,  the  plots  in  figure 
B-2  reveal  that  there  is  a  rather  broad  maximum  in  the  deflection,  meaning  that  the  choice  of 
ratio  r  is  not  critical. 

The  deflection  criterion  for  fourth-power  noise  is  depicted  in  figure  B-4.  For  low  input 
SNR,  the  best  ratio  r  is  approximately  0.65,  in  agreement  with  the  value  used  for  the  ROCs  in 
figures  21  through  23  in  the  main  report.  The  case  of  r  =  0.75  (used  in  figure  24)  has  virtually 
the  same  deflection  for  small  input  SNR  and  is,  in  fact,  slightly  better  than  r  =  0.65  as  the  SNR 
increases.  However,  the  broad  maxima  of  the  deflection  curves  indicate  that  the  exact  choice  of 
saturation  point  is  not  critical.  The  important  feature  of  the  nonlinearity  is  the  saturation  effect 
itself,  which  suppresses  the  heavy-tailed  character  of  the  non-Gaussian  noises  considered,  and 
allows  the  weak  input  signal  to  make  significant  contribution  to  the  output  over  many  samples  K. 


B-2 


%  deflection  ahn 


clear  all 
close  all  hidden 
Nu=200; 

Na=200; 

K=1000; 

Nb=ll; 
amax=4; 

q2=sqrt(2); 
q3=sqrt(3); 
qk=sqrt(K); 
p2=pi/2; 

a=[  1  :Na]/Na*amax; 
d=zeros(Na,Nb);  %  Deflection 

for  m=l:Na 
du=a(m)/Nu; 

u=[0:Nu]*du;  %  0  to  a(m)  Notice  limited  range. 

u2=u.^2; 

f=u2-a(m)^2;  %  Limiter  -  squarer 

%  pn=3/q2/pi./(l+2*u2)./(l+.5*u2);  %  (F) 

%  pn=  5*sech(p2*u);  %  (S) 

pn=sqrt(.5/pi)*exp(-.5*u2);  %  (G) 

yn=du*(2*sum(f.*pn)-f(  1  )*pn(  1 )); 
y2n=du*(2*sum(f.  *f.*pn)-f(  1  )^2*pn(  1 )); 
sn=sqrt(y2n-yn^2); 
forn=l:Nb 
dB=-10+(n-l)*.5; 

sigs=10^(dB/20);  %  Signal  standard  deviation 
w=q3*sigs;  %  Flat  signal  PDF 

%  cnp=.5+(2*atan((u+w)*q2)-atan((u+w)/q2))/pi;  %  (F) 

%  cnm=.5+(2*atan((u-w)*q2)-atan((u-w)/q2))/pi;  %  (F) 


%  cnp=2/pi*atan(exp(p2*(u+w)));  %  (S) 

%  cnm=2/pi*atan(exp(p2*(u-w)));  %  (S) 

cnp=phi(u+w);  %  (G) 

cnm=phi(u-w);  %  (G) 

psn=.  5/w*(cnp-cnm); 


ysn=du*(2*  sum(f.  *psn)-f(  1  )*psn(  1 )); 
d(m,n)=qk*(ysn-yn)/sn; 
end 
end 

figure(l),  elf 
plot(a,d,'k','linewidth',.7) 
grid  on,  zoom  on 
orient  tall 

Figure  B-1.  Program  for  Deflection  Criterion 


%  Number  of  u  increments 
%  Number  of  a  values 
%  Number  of  time  samples 
%  Number  of  SNRs 
%  Abscissa  range 


B-3 


2  2.5  3  3.5 

r  =  b/o 

Gaussian  Noise  and  Processor  (B-6) 
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